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Abstract 

The ^^Nb and "^^Sc NMR spectra in PbSci/aNbi/aOa (PSN) and 
PbMgi/3Nb2/303 (PMN) disordered relaxor ferroelectrics at the temperature 
T > To (To is the temperature of the dielectric susceptibihty maximum) 
have been studied. Spectra analysis was performed on the base both of the 
analytical description of NMR lines shapes, allowing for homogeneous and in- 
homogeneous broadening related to a random distribution of the electric field 
gradients and numerical Monte Carlo method taking into account electric field 
gradients originated from random distribution of Mg, Sc and Nb ions (which 
may be shifted or not) over B-type cation sites. 

The observed 1/2 < — > -1/2 transition spectrum both of the ^^Nb and 
^^Sc nuclei in the PSN was shown to contain a narrow (3-4 kHz) almost 
isotropic part and a broad strongly anisotropic part. These two components of 
NMR spectra are related to 1:1 Sc/Nb ordered and compositionally disordered 
regions of the crystal, respectively. It was shown that in the disordered regions 
Sc^^, Nb^''"and ions are shifted from their cubic lattice sites at one of 
three possible directions: < 100 >, < 110 > or < 111 >. 

In PMN the NMR spectrum of ^^Nb contains practically only the broad 
component. The portion of unbroadened spectrum that may correspond to 
ideal 1:2 regions accounts only for 1-2 percent of the total integral intensity. 
No evidance was obtained about existence of the 1:1 regions in PMN. The 
NMR data demonstrate that in PMN the cubic symmetry at T > Tq is locally 
broken due to ions shifts similar to that in disordered PSN. The values of the 
ion shifts were estimated in the point charges point dipoles approximation of 
the electric field gradients calculation both in the PSN and PMN. 
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I. INTRODUCTION 



Relaxor ferroelectrics |1|J^ have been attracting considerable attention in recent years 
due to their unusual physical behavior. Lead magnesium niobate, PbMgi/3Nb2/303 (PMN), 
is probably the best investigated perovskite relaxor crystal. Many of novel models which 
try to describe the physical properties of PMN are grounded on the assumption of the 
existence of a 1:1 Mg/Nb ordered microregions at least in 1/3 of the volume of a crystal. 
On the other hand, the latest x-ray and neutron diffraction measurements performed on the 
PbSci/2Nbi/203 (PSN) with different degree of cation ordering show that ferroelectric phase 
transition undergoes in both ordered (1:1 Sc/Nb) and completely disordered samples 0. It 
might mean that long-range ordering can not be suppressed completely by compositional 
(chemical) disorder. Moreover, a large shifts of all type of ions from cubic regular positions 
in PMN as well as in PSN already at T > Tq (Tq represents the temperature of dielectric 
constant maximum) were obtained from the computer treatment of the x-ray and neutron 
diffraction peaks. Due to this, the phase transition even in PSN becomes incomprehensible 
because, as it follows from x-ray data, at T > To, for instance, Pb^"*" ions are shifted away 
of about 0.027 nm with equal probability along one of the directions < 100 >, < 110 > 
or < 111 >. This consideration left unclear, however, if these ion shifts are the static or 
dynamical ones. 

Another difficulty in the interpretation of the x-ray and neutron diffraction data in relax- 
ors is related to amorphous character of observed diffraction peaks even at T << Tq. These 
methods provide an averaged view of the structure resulting in the cubic (PMN) or close 
to cubic symmetry (PSN). Broadening of diffraction peaks is caused by the compositional 
disorder of these materials. Thus, local methods of magnetic resonance and, in particular. 
Nuclear Magnetic Resonance (NMR) can be extremely useful in the study of such disordered 
systems because in NMR experiment the nuclei are sensitive to their environment only at 
a distance less then 1-2 nm. It provides also possibility to distinguish between the static 
and dynamic ion shifts that is practically unattainable in x-ray diffraction methods. 

In our previous papers (see, for instance, Ref. @]) we have already studied ^^Nb and 
^^Sc NMR spectra both in the PMN and PSN. In particular, it was shown that at T >> 
Tq NMR spectra of ^^Nb are mainly determined by the electric-field-gradient contribution 
originated from disorder in Nb/Mg(Sc) lattice sites. From NMR spectra analysis evidences 
were obtained about the existence of ordered and disordered microregions in both materials. 
Unfortunately, earlier measurements were carried out mainly at a low Larmor frequency 
[i>L ~ 49 MHz) so a bigger part of the ^^Nb spectrum belonging even to 1/2 < — > "1/2 
transition remained invisible and could not be analysed. More detailed ^^Nb and ^^Sc 
NMR spectra were recently presented in Ref. and 0. However, only the temperature 
dependence of the second moment of the narrow subset of NMR lines has been analysed. A 
broad background observed in the spectrum was attributed to unresolved satellite transitions 
and it was not analysed. 

The purpose of the present paper is to study local structure both of PMN and PSN at 
temperatures T > To on the base of ^^Nb and ^^Sc NMR measurements performed at Larmor 
frequency — 93 MHz. We believe that complete 1/2 < — > -1/2 transition spectrum can 
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be observed in such conditions. Because one of studied materials (PSN) belongs to a 1:1 
system, where the 1:1 Sc/Nb order really exists it should be thus interesting to see what 
are the common features and what are the differences in NMR spectra of these two types of 
relaxors. 

Description of the NMR line shape was performed on the base of analytical formulas, 
allowing for both homogenious and inhomogenious mechanisms of line broadening. As a 
second method of NMR spectra analysis we used numerical Monte Carlo simulation, which 
allowed to take into account random distribution of Mg, Nb and Sc ions (shifted away or 
not) over B-type sites of AB'B"03 perovskite lattice. 

The structural data obtained are compared with analogous data resulting from x-ray 
diffraction methods. Finally, the problem of 1:1 Mg/Nb ordered microregions in PMN is 
discussed as well. 

II. SAMPLES AND EXPERIMENTAL DETAILS 

The measurements were carried out on single-crystal samples of PMN and PSN grown 
on a seed crystal from solution in melt. Samples had a good optical quality. The sample 
dimensions of PSN were (4x4x4) mm^, and PMN - up to (8x8x8) mm^ with the surfaces 
parallel to crystallographic (001) planes. 

The NMR measurements were performed in a 9.2 T superconducting magnet at "^^Sc 
Larmor frequency pl = 92.3 MHz (I = 7/2) and ^^Nb ul = 92.9 MHz (I = 9/2). Tem- 
perature was stabilized in a continuous-flow cryostat with an accuracy of about 0.1 K. The 
inhomogeneously broadened spectra of ^'^Nb and ^^Sc (1/2 < — > -1/2 transition) have been 
obtained using the Fourie transform of the spin echo after a 90^: - r - 90^ - r - acquisi- 
tion pulse sequence. The width of the 90° pulse was 2.1 fis and 2.4 fis for ^^Nb and ^^Sc 
resonances, respectively. A time delay r was usually 40 fis for both of nuclei. 

III. EXPERIMENTAL RESULTS AND THEIR INTERPRETATION 

A. s^Nb and "^^Sc NMR spectra 

Measurements were performed at T = 400 - 450 K because we were interested mainly in 
the local structure of PMN and PSN in high temperature phase, i. e., at T > Tq (Tq ~ 355 
K and 310 K for PSN and PMN, respectively). As an example in Fig. |I] we presented ^^Sc 
NMR spectra measured for three characteristic crystal orientations in the external magnetic 
field: B II [001], [Oil] and [111] . The characteristic feature of the spectrum for B II [001] 
and [111] is its strong asymmetry while for B II [Oil] the spectrum is almost symmetric. 
Similar spectra of ^^Nb in the PSN and PMN are shown in Figs. |^ and ^ respectively. Let 
us discuss first the Sc and Nb spectra in the PSN. 

Both the Sc and Nb spectra have approximately similar structure and both contain a 
narrow component of the width of (3 - 4) kHz and a wide component that is different for 
the Sc and Nb nuclei. The narrow part of spectra is well described by a Gaussian and is 
completely isotropic. On the contrary, the wide part is essentially anisotropic and can not 
be precisely described by a simple function. In rough approximation, however, it can be also 
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described by a Gaussian (see in Ref. 0). In such approach the ratio of integral intensity of 
narrow and wide parts of the spectrum is approximately equal (40:60) and does not depend 
on the crystal orientation. Note that due to large width of the ^^Nb spectrum its wide 
component has very small intensity and thus its anisotropy manifests itself much less than 
that in the ^^Sc spectrum. However measurement of whole the ^^Nb spectrum in PSN, using 
the sweep of the irradiation frequency, indeed showed that the Nb spectrum shape is very 
close to that for ^^Sc (insert to Fig. 

^^Nb spectrum in PMN is somewhat complicated. In this spectrum we can conditionally 
separate three components: (1) low-intensity isotropic line with the width of 3 kHz; (2) 
spectrum around with the width of 15 - 18 kHz; (3) wide and strong anisotropic compo- 
nent like that in PSN but more intensive. The narrow part with Az/ = 3 kHz accounts only 
for about 1.5-2 percent of the total integral intensity of the spectrum. As can be seen from 
Fig. 1^, beside the narrowest component the whole spectrum is strongly anisotropic. 

The anisotropy of spectra as well as asymmetric form with respect to for both of 
^^Nb and ^^Sc nuclei and in both the PMN and PSN points out that given spectra belong to 
1/2 < — > -1/2 transition the frequency of which is shifted by the second order quadrupolar 
contributions |^. This quadrupole contribution exists only for nuclei with spin I > 1/2 and 
is related to the interaction of the nuclear quadrupole moment with the electric field gradient 
(EFG). Large values of quadrupole moments of ^^Nb and ^^Sc niclei {eQ = — 0.28e ■ 10~^® 
m^ and — 0.22e ■ 10~^^ m^ for ^^Nb and ^^Sc, respectively) usually result in predominance of 
the quadrupole contribution to NMR frequency shift with respect to other possible sources, 
as chemical shift and magnetic dipole-dipole interaction. Therefore we can draw conclusion 
that main features of ^^Nb and ^^Sc NMR spectra in both the PMN and PSN are related to 
characteristic features of their EFG, that depend on a lattice ions charges and their positions, 
i.e. local structure of crystals. 



B. Analysis of ^^Sc and ^^Nb spectra in PSN 

The integral intensity ratio (40:60) obtained above for the narrow and wide parts of the 
"^^Sc and ^^Nb NMR spectra in PSN is in agreement with the Sc/Nb degree of ordering which 
varied in the studed crystals in the range of 30 - 40 percent. So, it can be assumed that 
narrow and wide parts of the spectra belong to the Sc and Nb nuclei situated in the ordered 
1:1 and disordered regions of the crystal, respectively. Actually, in the ordered regions Sc 
and Nb ions are surrounded by a symmetrical ion configuration. Therefore in such sites 
the electric-field-gradients and anisotropic chemical shift must be close to zero and we can 
expect a NMR line broadened only by the magnetic dipole-dipole interaction. Calculation 
of this dipole-dipole width gives the value of about 3 kHz supporting the origin (ordered 1:1 
regions) of narrow lines in ^^Sc and ^^Nb NMR spectra. 

Evidently, the wide part of NMR spectra belongs to resonances from nuclei located in the 
disordered regions of the crystal, where random distribution of nonequal charges B'/B" ions 
produce electric-field-gradients. It follows from recent x-ray and neutron diffraction data [Q] 
that in these regions the lead and oxygen ions may be shifted away from their regular cubic 
positions of about 0.01 to 0.035 nm, which gives an additional contribution to the EFG 
tensor components. Of course, due to disorder in the ions positions we can expect broad 
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distribution in values of the EFG tensor components which manifestes in inhomogeneous 
broadening of NMR hues. 

The resonance frequency shift related to quadrupole mechanism for the considered case 
of ± 1/2 transition is proportional to the squared EFG. It can be written as follows 



^i% = §-^{ni+^)~l)m^), (1) 



where ur, = '^'^ ea = V n = '"'■'^''•'yy 
The function fn{6,ip) in Eq.(|lD describes the dependence of frequency shift on mutual 
orientations of external magnetic field and EFG tensor principal axes, vl is the Larmor 
frequency, x,y,z form the pricipal EFG tensor reference frame. 

Angular variation of the wide part of NMR spectra can be qualitatively understood if to 
assume presence of an axial EFG pointed at <001> cubic directions. In this case angular 
part of expression (0) takes a simple form: 

/(^, ^) = (1 - cos2 e) (9 cos2 e~i). (2) 

o 

Expected resonances for all six <001> directions of the EFG tensor are depicted as 
hystograms in Fig. ^ In this our simple interpretation we have to take into account that due 
to large disorder in the ions positions EFG values become randomly distributed that leads 
to strong inhomogeneous broadening of NMR lines. As a result, only the shift of the center 
of gravity is visible in the experimental spectrum under crystal rotation. Obviously, that 
because of nonlinear relation between resonance frequency shift and quadrupole frequancy 
inhomogeneously broadened NMR line acquires essentially asymmetrical form, which besides 
depends on orientation of a crystal in external magnetic field. For such complex spectrum its 
decomposition on simple Gaussian components and the analysis of their angular dependence 
can not give true values of quadrupole frequencies. 

In what follows, we introduce the distribution of largest EFG component which leads to 
the distribution of quadrupole frequencies pq. Let us suppose that this distribution function 
can be expressed in Gaussian form, i. e. 



where Vq is the mean value of the quadrupole frequency and A is the width of its 
distribution. 

Allowing for the nonlinear relation between u and vq (see Eq.(|l|)), namely 



^ = «4, (4) 



the shape of inhomogeneously broadened line can be written as (Ref. 
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Substitution of Eq. into Eq. gives 



,0 \2 



exp( 



2aA2 



cosh — . 
A2 V a 



(6) 



To take into account the additional homogeneous broadening mechanism with Gaussian 
hneshape (dipole-dipole interactions) one has to make the integration 



1 f'^ 
/(^) = / exp 
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(7) 



where b being a half of the Gaussian width. 

One can see that Eqs. (0-0) give an asymmetrical lineshape which depends also on the 
parameter a(^,(y9) that takes into account the mutual orientations of EFG and magnetic 
field Bq. 

One can expect that in disordered materials there are microregions with different "degree 
of order" and different directions (and mean values) of the ion displacements. Obviously, 
in such a case the NMR line should be the composition of the lines "stemming" from these 
microregions, i.e. 



(8) 



where subscript "i" numerates the microregions. 

The results of simulation on the base of Eqs.(^^ of the "^^Sc NMR line shape for B|| [001] 
are presented by solid line in Fig. |^(a) and (b). Figure ^(a) shows the line shape for the 
case of EFG with largest principal axis points in (001) directions only. The contribution 
from ordered regions, where i/q = was included also. The best fit has been achieved for 
I t'g I = 755 kHz and A = 210 kHz. The contribution of homogeneous broadening was 
estimated as 5 = 2 kHz. One can see that calculated NMR line has a strongly asymmetrical 
shape, which fits very good the right "shoulder" of the measured spectrum. The fitting of 
left "shoulder" is not so good. This can be expected because for (001) orientation of EFG 
only a positive quadrupole shift of resonance frequency exists (see Eq.(4) when d = 90°; 0°). 
Therefore, to fit the left shoulder of the observed spectrum the other orientations of EFG 
axes are required. The negative shift is possible for EFG axes pointing along (Oil) and 
(111) directions. The reason for this type local distortions can be related to substitutional 
disorder in the relaxors. 

Really, the probability of configurations with k ions of B type in the nearest neighbour- 
hoods of B or B ions can be written as 



6!i/^(l -y) 
k\{6-k)\ 



6-fc 



2/ = 1/3 



(9) 



One can see that P^ for k 



1, 2, 3 are close to each other, while the probability of having 
A; = 6 or is much smaller than the most probable value P^ = 0.328 corresponding to k 
= 2. It is easy to see that configurations with A; = 1,2,3 can be the response of the three 
afore-mentioned lattice distortions. 
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Finally, the experimental line shape of Sc was well fitted by the following theoretical 
curve 

Ig (U) = Colord (l^) + C^l/<001> {f^) + C2/(110> {f^) + C3/<111> (10) 

with Co ^ 0.3, Ci ~ 0.45, C2 ~ 0.11, C3 ~ 0.14; the subscripts denote the direction of 
the EFG axes. The result of ^^Sc complete spectrum simulation is presented in Fig. ^(b). 

The obtained set of fitting parameters {uq and A are presented in Table ID permits to 
describe theoretically the NMR line shapes measured for any other orientation of magnetic 
field. In particular in Fig. ^(c) and (d) the observed and calculated line shapes for B||[011] 
and B||[lll] are shown. One can see the good agreement between measured and calculated 
line shapes. 

Similar simulations performed for ^^Nb NMR spectra are presented in Fig. ^ and cor- 
responding fitting parameters are listed in Table |. The agreement between measured and 
calculated spectra seems also good, however, due to small intensity of the wide part of the 
spectrum obtained parameters may contain significantly bigger miscalculation than that for 
the "^^Sc spectrum. 

C. Analysis of ^^Nb NMR spectra in PMN 

It can be seen from Fig. ^ that practically all the ^^Nb spectrum in PMN is strongly 
broadened in comparison with the dipole-dipole width that is equal to 2.8 - 3.6 kHz (as 
in the case of PSN). The portion of the spectrum unbroadened by quadrupole mechanism 
or by anisotropic chemical shift accounts only for 1-2 percent of the total intensity. It is 
clear that the Nb nuclei, which give the contribution into this narrow spectrum, locate at 
the positions close to the regular cubic ones. Possible origin of these crystal regions will be 
discussed in the next section. 

The wide, main part of the ^^Nb spectrum is related to the disordered, locally noncubic, 
regions of the crystal , as in the case of PSN. Contrary to PSN, the portion of such regions 
in PMN accounts for 98% of the total volume of the crystal. 

^^Nb spectrum for disordered PMN was fitted assuming local distortions along <001>, 
<110> or <111> directions, i.e. as in the case of the PSN by Eq. (||). Because the 
procedure of the simulation is completely similar to that described for PSN we confined here 
only presentation in Table | of spectral parameters derived. Compareing spectral data for 
disordered PSN and PMN one can emphesize a similarity in the symmetry and sometimes 
even in the values of electric field gradients at Nb positions. One can see also that dispertion 
in PMN is larger than in PSN. This is in agreement with the fact that in PMN local lattice 
distortions should be much larger due to bigger difference in the ionic charges of B' and B" 
cations. 

Summarizing, three main conclusions should be made. 

(i) In the ordered regions, which account up to 30% in the PSN and 2-3% in the PMN 
of the total crystal volume, lattice ions occupy cubic lattice positions where electric field 
gradients are close to zero. 

(ii) In the disordered regions there are non-zero EFG along pure cubic directions (<001>, 
<110> or <111>). Because both the PSN and PMN has in average cubic symmetry at T > 
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To, these non-zero EFG can be related to randomly distributed unequal charges of B-type 
cations and/or random shifts of the ions from their positions in the cubic lattice. In this 
case the essential part of ions (~ 30-55%) experiences shift along <001> axes. 

(iii) EFG tensor components are significantly distributed around their mean values that 
manifests in broad asymetric NMR line shape. 



In the previous section on the basis of analytical calculation of the shape of inhomoge- 
neously broadened NMR line the mean values of quadrupole coupling constants in Sc and 
Nb cation positions as well as their dispertions were obtained. The method applied by us 
has allowed to take into account only distribution of the quadrupole frequencies (i^q), but 
the presence of random deviations of EFG axes from their most probable directions, i.e the 
possible fluctuations of fniO.ip) (see Eq. (|I|)) remained discounted. Note that this effect is 
mainly related to a distant ion contributions. Ignoring of these fluctuations has resulted in 
too much high values of homogeneous width, which was accordingly 4 kHz and 8 kHz for 
the ^^Sc and ^^Nb in PSN and 18 kHz for the ^^Nb in PMN. Let's remind, that dipole width 
for both the ^^Sc and ^^Nb accounts only 2.8-3.2 kHz. 

The second reason on which we have undertaken numerical calculation of EFG is related 
to attempt to define ions displacements, which could be responsible for symmetry and mag- 
nitude of an EFG in both relaxor systems. Taking into account complexity and substantially 
uncertainty of local structure of these compaunds we used the simple model of point charge 
point dipole, which, however, took into account effects related to partial covalency of bonds 
in the BOg octahedron as well as electronic polarizibility of ions. Obviously due to limi- 
tations of the point charge model the obtained ions shifts should be considered only as an 
estimated ones, the degree of which rehabihty depends on correctness (luck) of a choice of 
such empirical parameters of ions as their effective charges and polarizibility. 

In the point charge point dipole approximation the electric field gradient components at 
the site of the nucleus are given by 



where and Hk are the electric point chages and dipole moment on the k-th ion. is 
the distance between the k-th ion and the observation point, and 7oo is the anti-shielding 
parameter 

For the conviencity three contributions into EFG were detached in the whole sum (|ll]). 
The first contribution iYip) was related to random distribution of charges of B'/B" ions at 
lattice points of ideal ABO3 perovskite structure. Contributions of ions at a distance up 
to five lattice constants (2 nm) were summed up. However, usually the distance up to 1.2 
- 1.5 nm was sufficient due to the rapid decay of Vj^ with the distance (see Eq. |TTp. In 
the calculations the B type ions were randomly distributed at lattice sites with the help of 
random number generator or they were located in the 1:1 or 1:2 ordered regions. The main 
condition for any ion configuration was electroneutrality, i. e. the stoichiometric (in average 
1:1 for PSN and 1:2 for PMN) composition in the calculated region was conserved. 



D. Numerical calculation of the EFG tensor and NMR line shape 




(11) 
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As an example, the result of the calculation of largest component Vzz as a function of 
polar angle 6 for the PMN is presented in Fig. |^. Here we took effective valence charges of 
Mg^"*" and Nb^"*" as one haft of their formal charges. It is seen that Vzz hes close or along 
cubic directions, but it is strongly distributed in both the magnitude and orientation. The 
mean value | Vzz \= 0.78*10^^ V/cm^ — ^^MHz) is much less than that expected from 
coupling constant measured. Thus, this contribution to the EFG can not be responsible for 
wide component of the ^^Nb NMR line. It can give only some contribution to homogenious 
broadening of NMR line because it leads to distribution of axes directions of the total EFG. 
Note that in PbSci/2Nbi/203 this part of EFG is neglectly small in comparison even with 
the dipole-dipole width, that is in good agreement with smaller difference between charges 
of Sc^^ and Nb^"*" and, in addition, 1:1 type relaxors have more symmetrical distribution of 
B'/B" ions. 

The second contribution (V^y ) into EFG took account for the shifts of B'/B" ions them- 
selves relatively to the nearest six oxygen ions. Possible shifts of Pb^"*" ions were not taken 
into account due to their smaller (10-15 times) contribution into EFG for the same values 
of the shifts. 

In a general case for arbitrary B type ion shift with respect to the center of the undistorted 
oxygen octahedron, the following expressions for EFG tensor components were obtained : 

(2) 42^ 

V^f = '-^{2dl-dl-dl) (12) 
21 ■ z^(^^ 



y(2) 

yy ^5 

y(2) 



21 ■ z"^^ 





-dl 


-< 




-dl 


-4 


{2dl 


-dl 


-< 
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where (d^,, dy, d^) are vector components of the Nb or Sc shift, Zq"' is the effective valence 
charge of the oxygen, 6 = a/2 (a is the lattice constant). For the simplisity the 2;™' was 
taken as an avarage value over all six oxygen ions. Its value may lay in the range between 
-1.63 I e I (BaTiOg ^) and -0.85 | e | (KNbOg Q). 

The last term in Eq. (|l^) is related to the induced electric dipole moments /Xj. Assuming 
that only electronic polarization of the nearest six oxygens gives marked contribution to the 
EFG (the other ions in the unit cell and the next nearest oxygens give contributions that 
are an order of magnitude smaller), the dipole contibution is: 

V„ = K,, = -^d (13) 

where /i is assumed to be in the direction of the Z axis of the field gradient. 
In the case of <110> and <111> symmetry of lattice distortions nondiagonal components 
of the EFG tensor have to be also appeared. They are given by expression 

^^. = -^^.' '^3- (14) 
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Electronic dipole moment dependes on ionic polarizability and value and direction of 
ions shifts and it can be usually calculated using knowledge of the polarization and internal 
electric field. Obviously for disordered relaxors such a way is not suitable because here is only 
local polarization and large uncertainty in the ions positions. Therefore, for the estimation 
of Vij we applied a method of effective dipole moments. In particular, it is known, that 
in perovskite ferroelectrics the spontaneous polarization can be well described by the Born 
effective charges: 



e 



(15) 



where Q is the cell volume, and the Born effective charges Z* are 0.82(K), 9.13(Nb), 



-6.58(Oi), -1.68(02) in tetragonal phase of KNbOg ^ and 3.87(Pb), 7.07(Ti), -5.71(Oi), 
-2.51(02) in PbTiOa [1^. The high values Z* indicate on larger contribution of electronic 
dipole part to whole polarization in the comparison with rigid ionic part, that should be as 
a rule in high polarizable perovskite ferroelectrics [0 . 

Taking into account a special role of Pb ions, namely their strong hybridization with 
oxygen ions, we used in our calculations the effective charges derived in PbTiOa. For sim- 
plisity the same avaraged value of Zq = —4 | e | was taken for evaluation of the electronic 
dipole moment fii = Zq * di in expressions (^^-|T^). One can easy see, that under these 



conditions expressions ([T^ ) and ([T3| - p!4D become similar, thus in the calculation of the whole 
dipole contribution to the EFG one can use only its electronic part, because Zq >> Z^K 

m 



and BaTiO'^ [16 



Note that similar results were obtained in KNbOa 

Taking into account that the ion shifts can have random values due to the disorder of the 
Mg, Nb and Sc positions, the ion shift distribution function was assumed in the Gaussian 
form: 



fid) 



1 



a(27r)V2 



exp 



jd - dp) 
2^2 



(16) 



where do is mean value of the shift and a is its dispersion. 

Complete NMR line shape was calculated taking into account electric-field-gradient con- 
tribution, magnetic dipole-dipole interaction, as well as chemical shift mechanisms. Since 
the contribution of the last mechanism was much less with respect to the quadrupole one 
it manifests itself only in the central narrow part of the spectrum. Its contribution was 
accounted by the renormalization of dipole-dipole linewidth. NMR spectrum intensity /(z/) 
was obtained by summarizing the contributions of the NMR lines from a large number of 
clusters in which the B type ions could be shifted in one of three directions: <001>, <011>, 
<111> with their resonance frequencies given by Eq. (||). 



N 

E 

n=l 



exp 



In 2) 
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(17) 



Here, Un is the resonance frequency for nucleus in the center of n-th cluster, 5 is a half of 
the dipole-dipole width and is the total number of clusters (usually, N = 15000 - 20000). 

In the calculation of second-order quadrupole shifts for ^^Nb the parameter = — 16 
was taken as in other ABO3 crystals like KNbOa and LiNbOa |lT5| , p!7[| . For the '^^Sc ion the 
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antishielding parameter is rather uncertain. However, clearly its value has to be smaller 
of about two times in comparison with ^^Nb due to smaller number of electronic orbitals. 
Value = —9 leads to the best agreement with the experimental spectrum. 

As an example, the results of ^^Sc complete spectrum calculation in PSN for several 
crystal orientations are presented in Fig. |^ and corresponding values of ion shifts as well as 
their dispersion in Table It is seen excelent agreement between measured and calculated 
spectra including frequency region around ul, where analytical method gave much worth fit. 

The results of calculation in PMN for two crystal-orientations (B II [001] and [Oil]) are 
presented in Fig. ^. The individual contributions which were related to the Nb shifts in 
different directions with respect to oxygen positions are presented as well. The values of ion 
shifts and their dispersion obtained by fitting of the experimental NMR spectra are listed in 
the Table |IT1[ Excellent similarity of ion shifts in PMN and disordered PSN is obtained. On 
the other hand, a large relative dispersion of niobium (or rather oxygen) shift is noted along 
<111> and <011> directions. Such a large fluctuations in ion shifts can show their dynamic 
character. Actually, we can not expect full motional narrowing of NMR spectrum, since the 
frequency shift of the -1/2 < — > 1/2 transition depends on the square of EFG. However, 
this assumption can be supported by measuring and analysis of NMR spectra changes with 
decreasing temperature. We plan to carry out such analysis in the nearest future. On the 
other hand, large dispersion of the quantities is characteristic feature of relaxors. 

Three main comments should be made. 

(i) The values presented in Table P-|nT|are the relative shifts between Nb(Sc) and oxygen 
ions, i. e., dm{Sc)-o- Larger contribution into dm{Sc)-o values due to oxygen ions shift rather 
than that of Nb(Sc) ions is expected in accordance with x-ray data [^. 

(ii) The real values of ion shifts can be somewhat different from those given in the 
Tables [rT[]ITT| due to simplifications allowed in the calculations, including point charges model. 
However the conclusion about existence of the displacements, their directions and dispersion 
is not related to the model of EFG calculation because it follows directly from observed shape 
and width of ^^Sc and ^^Nb NMR lines. 

(iii) In the ordered regions the Sc^"*", Nb^+and 0^~ ions occupy completely cubic lattice 
positions, i. e. they are not shifted. On the contrary, in the disordered regions these all ions 
are randomly shifted away from their cubic lattice sites in one of the directions: <001>, 
<011>or <111>. The larger part of ions (~50%) experiences shift along <001>axes. 



IV. DISCUSSION 



Evidently, basic question in the understanding of the origin of unusual physical properties 
of relaxor ferroelectrics is their local structure, because at any temperature and even on a 
micrometric scale the average symmetry is cubic rather than any lower one (see, for instance, 
I8| , p!9[l and references therein). Note, that the existence of diffusive scattering in the powder 



neutron diffraction patterns observed in PMN below 340 K was interpreted in Ref. ||2^ as 
the occurrence of local polar order inside the average paraelectric matrix. However, it is still 
a matter of debate whether this polar order is organized as small ferroelectric domains or as 
a dipole glass. In contrast to PMN, PSN at T < Tc (Tc ~ 350-360 K) undergoes into the 
i?3m rhombohedral ferroelectric phase PJ2T|. Moreover, this crystal can be easily obtained 



in 1:1 Sc/Nb ordered state. It is important to emphasise that the polar long range order 
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appears not only in perfect 1:1 Sc/Nb ordered PSN, but even in compositionally disordered 
crystal f^, i. e. the compositional disorder does not prevent establishing at least partly 
of a long-range order . To our mind this may correspond to mixed ferro-glass phase. It 
is important to compare PSN and PMN local structure on the base of NMR data because 
many of novel models, which explain the behaviour of PMN at low temperatures, are based 
on the existence of 1:1 Mg/Nb ordered micro-regions similar to 1:1 Sc/Nb regions in PSN. 

It is clearly seen from Figs. |I] and ^ that for 1:1 ordered regions in the PSN distinctive 
narrow line in both ^^Sc and ^^Nb NMR spectra are obtained. In these regions Sc and 
Nb ions are located in regular cubic lattice sites, where the total electric-field-gradient is 
zero or close to zero. The ^^Nb NMR spectrum in PMN also exhibits narrow component 
(Fig. 1^) but its integral intensity is only of about 1.5-2 percent of the total spectrum 
intensity (Tables [n|jlll|). Evidently, it is reason to ascribe this ^'^Nb resonance to the regular 
1:2 Mg/Nb ordered regions of the crystal. Note that this resonance could not belong to 1:1 
Mg/Nb ordered regions because in agreement with high resolution electron microscopy study 
the portion of such ordered regions (if they exists) should be approximately 30%. Thus, 
if local (~ 1.5 - 2.5 nm) 1:1 ordered clusters exist in PMN, the disagreement between NMR 
and electron microscopy data can be overcame in supposition that in these local regions 
niobium and oxygen ions are shifted as in completely disordered regions. This fact seems to 
be strange. Let us remind that in 1:1 Sc/Nb ordered PSN the Sc and Nb ions are situated 
in correct cubic lattice sites. Taking into account this argument, the model grounded on 
the alternation of Nb and (Mg2/3Nbi/3) layers which also explains the observation of 1/2 
<111> super lattice diffraction peak seems to be more convincing. In this model 

local strains induced by the charge disbalance between the Nb and Mg/Nb layers can be 
locally compensated by oxygen ion displacements which will give the contribution in the 
electric field gradient. Additional contribution to EFG is also due to (Mg2/3Nbi/3) layer 
compositional disorder. 

Obtained evidence about three possible directions of Nb ions shift relatively the sur- 
rounding ions (so that oxygen and lead ions shift can contribute also) seems to be very 
important. Physical nature of this phenomenon can be related to substitutional disorder in 
the relaxors A(B-^y3B2/3)03. 

We can not exclude also that ion shift in one out of three observed directions may be 
the result of ferroelectric phase transition which could be at Burns temperature T^ = 640 K 



25| . However with the temperature decrease strong random fields completely destroy this 
ferroelectric long range order because of correlation radius decrease at T < T^ (here 
is the correlation radius of the lattice that should be responsible for the phase transition 
at T^) As a result only a small movable polar clusters appear below T^. In such a 

picture strong dispersion and V-F law, observed in the temperature dependence of dielectric 



permitivity in PMN (see, e.g. ||2^ and ref. therein) can be considered as manifestation 
of reentrant phase (like that in magnetic systems |^) of aforementioned Burns reference 
phase. 

To clear up all aforementioned questions precise measurements of NMR spectra in the 
vicinity of the Burns temperature are extremely desirable. 

Finally, we would like to point out that proposed methods of NMR line shape analysis 
can be applied to other disordered materials where the NMR frequency shift is a nonlinear 
function of the random electric field gradients. 
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TABLES 



TABLE I. ^Sc and ^^Nb quadrupole coupling constants and their dispertion measured in 



disordered PSN and PMN at T = 450 K. 



Compound 


^^^'^ (MHz) 


A (MHz) 


1 1 (lO^Oym-^) 


Orientation 
of Yzz 


PbMgi /QNbo/QO'? 










45 So 


10.6 


2.9 


20 


<001> 




Q 1 




1 7 






8.8 


4.9 


16 


<111> 


93Nb 


46 


15 


68 


<001> 




42 


20 


62 


<011> 




56 


27 


83 


<111> 


PbMgi/3Nb2/303 










93Nb 


46 


24 


68 


<001> 




60 


29 


89 


<011> 




36 


29 


53 


<111> 



lABLE II. PMrame(x;rs of ions t^liifls in PSN al T = 420 K derived from NMR data. 





1:1 ordered 


<001> 


direction of ions shifts 
<011> 


<111> 


portion (%) 


30 


35 


15 


18 


dsc-o (nm) 





0.015 


0.016 


0.017 


a (nm) 





0.001 


0.001 


0.001 


dm-O (nm) 





0.016 


0.018 


0.022 


a (nm) 





0.005 


0.005 


0.005 



TABLE HL Parameters of ions shifts in PMN at T = 450 K derived from NMR data. 





1:2 ordered 


<001> 


direction of ions shifts 

<011> 


<111> 


portion (%) 


2 


56 


21 


21 


dNb-O (nm) 





0.016 


0.021 


0.017 


(7 (nm) 





0.005 


0.01 


0.01 



14 



FIGURES 



FIG. 1. Sc NMR spectra in PSN. The diagrams schematecally show expected NMR reso- 
nances for second order quadrupole contribution with EFG axis pointed at < 001 > directions. 

FIG. 2. ^^Nb NMR spectra in PSN. The inset shows the ^^Nb spectrum measured by sweep of 
the irradiation frequency. 

FIG. 3. '^^Nb NMR spectra in PMN. 

FIG. 4. Comparison between experimental (points) and calculated (solid line) ^^Sc NMR line 
shapes in PSN for a) and b) B||[001]; c) B||[011]; d) B||[lll]. Calculated line shapes include 
contribution from ideal structure regions plus axial EFG pointed at a) (001) directions; b), c) and 
d) (001), (Oil), (111) directions. 

FIG. 5. Comparison between experimental (points) and calculated (solid line) ^^Nb NMR line 
shapes in PSN for a) B||[001]; b) B||[011]. 

FIG. 6. Polar plot of the | V^z \ EFG component at Nb sites in PMN calculated for Mg^+ and 
Nb^"*" random distribution at lattice points of ideal ABO3 structure. 

FIG. 7. Measured and calculated ^^Sc NMR spectra in PSN at orientation (a) B II [001] and 

(b) B II [Oil]. Separated contributions into NMR spectrum connected with different type of ion 
shifts are shown in the bottom panels by: solid line (dsc-O II <001>), dash line (dsc-O H <011>) 
and dot line (dsc-o H <111>); contribution from 1:1 Sc/Nb ordered regions is shown by dash dot 
line. 

FIG. 8. Measured and calculated ^^Nb NMR spectra in PMN at orientation (a) B II [001] and 
(b) B II [Oil]. Separated contributions into NMR spectrum connected with different type of ion 
shifts are shown in the bottom panels by: solid line (d]V6-o H <001>), dash line ((ijv6-o H <011>) 
and dot line (djvb-o H <111>); contribution from 1:2 Mg/Nb ordered regions is shown by dash 
dot line. 
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figure 2 
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figure 3 
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